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Abstract. We analyze the influence of the magnetic field in the convexity properties 
of the relativistic magnetohydrodynamics system of equations. To this purpose we 
use the approach of Lax, based on the analysis of the linearly degenerate/genuinely 
non-linear nature of the characteristic fields. Degenerate and non-degenerate states 
are discussed separately and the non-relativistic, unmagnetized limits are properly 
recovered. The characteristic fields corresponding to the material and Alfven waves 
are linearly degenerate and, then, not affected by the convexity issue. The analysis 
of the characteristic fields associated with the magnetosonic waves reveals, however, a 
dependence of the convexity condition on the magnetic field. 

The result is expressed in the form of a generalized fundamental derivative written 
as the sum of two terms. The first one is the generalized fundamental derivative in the 
case of purely hydrodynamical (relativistic) flow. The second one contains the effects 
of the magnetic field. The analysis of this term shows that it is always positive leading 
to the remarkable result that the presence of a magnetic field in the fluid reduces the 
domain of thermodynamical states for which the EOS is non-convex. 
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There are many astrophysical scenarios governed by relativistic magnetohydrodynamical 
processes as, e.g., the production of relativistic jets emanating from Active Galactic 
Nuclei, the structure and dynamics of pulsar wind nebulae, the mechanisms triggering 
the explosion in core-collapse supernovae, or the production of Gamma Ray Bursts. 
These scenarios are nowadays the subject of intensive research by means of numerical 
simulations thanks to recent advances in numerical relativistic magnetohydrodynamics 
(RMHD) that exploit the fact that the RMHD equations obeying a causal equation of 
state (EOS) form a hyperbolic system of conservation laws pj. 

Matter at densities higher than nuclear matter density can undergo first-order phase 
transitions to various phases of matter, such as pion condensates [2], hyperonic- matter 
[3] or deconfined quark matter [U [5] . Several authors 0 0IE] have studied, from 
different points of view, the influence that those exotic states of matter at extreme high 
densities have on, e.g., the dynamics of stellar core collapse supernovae, the evolution 
of proto-neutron stars, or the collapse to black hole. 

The classical Van der Waals (VdW) EOS is a well known example of EOS displaying 
a first-order phase transition. Fluids having a thermodynamics governed by a VdW- 
like EOS exhibit, outside the region of the phase transition, non-classical gasdynamic 
behaviours in a range of thermodynamic conditions characterized by the negative value 
of the so-called fundamental derivative, Q ® EDI HU 

d 2 p 


G: = 



dV 2 
dp 


( 1 ) 


dV 


S 


p being the pressure, V := 1/p the specific volume (p is the rest-mass density) and s the 
specific entropy. The fundamental derivative measures the convexity of the isentropes 
in the p — V plane and if Q > 0 then the isentropes in the p — V plane are convex, 
leading to expansive rarefaction waves (and compressive shocks) [12]. In a VdW-like 
EOS, or in general in a lion-convex EOS, rarefaction waves can change to compressive 
and shock waves to expansive depending on the specific thermodynamical state of the 
system. These non-classical phenomena have been observed experimentally and their 
study is, currently, of interest in many engineering applications mm- 

Besides this thermodynamical interpretation of convexity, there is an equivalent 
definition due to Lax [15] that connects with the mathematical properties of the 
hyperbolic system. According to Lax’s approach, a hyperbolic system of conservation 
laws 0 is convex if all its characteristic fields are either genuinely non-linear or linearly 


1 The books by LeVeque m and Toro HZ] are recommendable references for those readers interested 
in the basic theory of hyperbolic systems of conservation laws. The monograph of [18j on finite- 
volume methods for hyperbolic problems pays special attention to non-convex flux functions (see their 
Sects. 13.8.4 -definitions of genuine non-linearity and linear degeneracy, and their relationsphip with 
convexity-, and 16.1 -devoted entirely to the study of scalar conservation laws with non-convex flux 
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degenerate. A characteristic field A is said to be genuinely non-linear or linearly 
degenerate if, respectively, 


CD 

Ik 

> 

(2) 

V U A • r = 0, 

(3) 


for all u, where V U A is the gradient of A(u) in the space of conserved variables, r is 
the corresponding eigenvector, and the dot stands for the inner product in the space of 
physical states. 

In a non-convex system, non-convexity is associated with those states u for which 
the factor V corresponding to a genuinely non-linear field, Eq. ([2]), is zero and changes 
sign in a neighbourhood of u. 

A virtue of Lax’s approach is that it can be applied to other hyperbolic 
systems in which the convex or non-convex character of the dynamics is governed 
by other ingredients beyond the EOS. Among these systems are those of relativistic 
hydrodynamics (RHD) and classical magnetohydrodynamics (MHD). In these two cases, 
the convexity of the system has been characterized with the sign of a generalized 
fundamental derivative that includes an extra term depending of the local speed of 
sound (in the case of RHD [19]) and the magnetic field (in the case of MHD [20]). 

In this work we use the approach of Lax to characterize, from a theoretical point of 
view, the effects of magnetic fields in the convexity properties of the RMHD system of 
equations as a previous step to explore its possible impact in the dynamical evolution 
of different astrophysical scenarios. The result is presented in the form of an extended 
fundamental derivative whose sign determines the convex/non-convex character of the 
RMHD system at a given state. Our result recovers the proper non-relativistic and 
unmagnetized limits. 

The paper is organized as follows. In Sect. 2, the equations of RMHD are introduced 
as a hyperbolic system of conservation laws. The transformation between primitive and 
conserved variables are explicitly written. In Sect. 3 the characteristic structure of the 
RMHD equations is discussed and the analysis of convexity in non-degenerate states 
presented. In Sect. 4 the analysis of convexity is extended to degenerate states. The 
non-relativistic, unmagnetized limits are recovered in Sect. 5. Section 6 includes a short 
summary and presents the conclusions. Finally, there is an Appendix that displays 
the Jacobian matrices of the RMHD system in quasi-lincar form, necessary for the 
characteristic analysis of Sect. 3. 

2. The equations of ideal relativistic magnetohydrodynamics 

Let J M , and be the components of the rest-mass current density, the energy- 

momentum tensor and the Maxwell tensor of an ideal (infinite conductivity) magneto¬ 
functions-). 

§ Throughout this paper, Greek indices will run from 0 to 3, while Roman run from 1 to 3, or, 
respectively, from t to z and from x to z 1 in Cartesian coordinates. 
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r = pvR (4) 

T^ v = ph*u»u v + g lw p* - W (5) 

*F ftv = uW - u v b **, ( 6 ) 

where p is the proper rest-mass density, h* = 1 + e + p/p + b 2 / p is the specific enthalpy 
including the contribution from the magnetic held ( b 2 stands for b^b^), e is the specific 
internal energy, p is the thermal pressure, p* = p + b 2 / 2 is the total pressure, and g 
is the metric of the space-time where the fluid evolves. Throughout the paper we use 
units in which the speed of light is c = 1 and the (47T) 1 / 2 factor is absorbed in the 
definition of the magnetic held. The four-vectors representing the fluid velocity, tC, and 
the magnetic held measured in the comoving frame, b satisfy the conditions = — 1 
and u^b^ = 0. 

The equations of ideal RMHD correspond to the conservation of rest-mass and 
energy-momentum, and the Maxwell equations. In a hat space-time and Cartesian 
coordinates, these equations read: 


= 0 

(7) 

T ^ = 0 

(8) 

= 0, 

(9) 


where subscript ( )/t ) denotes partial derivative with respect to the corresponding 
coordinate, (t,x,y,z), and the standard Einstein sum convention is assumed. 

The above system can be written as a system of conservation laws as follows 

<9U dF‘ 

dt dx l 

where V = (p, v l ,e, B l ) T is the set of primitive variables. The state vector (the set of 
conserved variables) U and the fluxes, F'', are, respectively: 



( D\ 


\ Bl ) 

( D v l 

S j v i + p*8 ij - b'lf/W 
rv i + pV - b^B^W 

v ip>k _ ykj^l 


( 11 ) 


( 12 ) 
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In the preceding equations, D , Si and r stand, respectively, for the rest-mass 
density, the momentum density of the magnetized fluid in the j-direction, and its total 
energy density, all of them measured in the laboratory (i.e., Eulerian) frame: 


D = pW, 

(13) 

S l = ph*W 2 v i - b°b\ 

(14) 

T = ph*W 2 -p* - {b 0 ) 2 - D. 

(15) 


The components of the fluid velocity trivector, v l , as measured in the laboratory frame, 
are related with the components of the fluid four-velocity according to the following 
expression: u M = W(l,v l ), where W is the flow Lorentz factor, W 2 = 1/(1 — v l Vi). 

The components of the magnetic field four-vector in the comoving frame and the 
three vector components B l measured in the laboratory frame satisfy the relations: 


b° = Wv k B k , 

V = U + b°v\ 

W 

Finally, the square of the modulus of the magnetic field can be written as 


b 2 = 


B k B k 

W 2 


+ (' v k B 


k \2 


(16) 

(17) 

(18) 


The preceding system must be complemented with the time component of 
equation (j2j), that becomes the usual divergence constraint 

dB i 


dx i 


= 0 . 


(19) 


CLg . 


An EOS p = p(p, e) closes the system. Accordingly, the (relativistic) sound speed 
, e being the mass-energy density of the fluid e = p( 1 + e), satisfies 

dp 


dp 

de 


p dp 

ha 2 = x + — «, with y := — 
p z dp 


and k := 


de 


3. Characteristic structure of the RMHD equations and analysis of 
convexity in non-degenerate states 

The characteristic information of the system of RMHD (fTOl) is contained in the set of 

<9F* 

eigenvalues and right eigenvectors {A a , r Q }® =1 of ( k B k , where B l := are the Jacobian 
matrices of the vectors of fluxes along the coordinate directions, and Q is an arbitrary 
unitary 3-vector. 

Since the dependence on U of the fluxes F* is implicit, it is useful to write the 
Jacobian matrices B l in terms of matrices involving only explicit derivatives with respect 
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to the primitive variables, V. If we define A 0 := ——, and A 1 := ——, then we have 

o V' o\ ' 

that B‘ = A 1 {A 0 )" 1 . Now, the sets of eigenvalues and right eigenvectors of the system 
in conservation form, {A a ,r Q }® =1 , and of the system in quasi-lincar form, {A*,r*}® =1 , 
satisfying (CkA k — A*M°)r* = 0, are related according to {A a ,r Q }® =1 = {A*,M°r*}. 
Matrices A 0 and (kA k are displayed in the Appendix. 

Once the eigenvalues and eigenvectors are known, we can analyze the convexity of 
the system studying the expression V a = VuA Q • r Q (see the Introduction). Finally, we 
can take advantage of the fact that, since A 0 is non-singular, then V a 7 ^ 0 if, and only 
if, Vf := VvA a • r* 7 ^ 0 , and perform the analysis of convexity in terms of Vf. 

The eigenvalues A Q are the solutions of the following polynomial expression for A 

Aa (Sa 2 - £ 2 ) (( 6 2 + pha 2 s )a 2 G - W s ~ 2 pha 4 - a 2 G£ 2 ) = 0, (20) 

where £ := ph + b 2 , Wf 2 := 1 — a 2 and quantities a, G and B were defined in ref. [T], 
a := 4> a u a , G := 4> a (j) a , B := (j) a b a , being, in our case, (f) a := (— \,(i) the normal to the 
wavefront propagating with speed A in the spatial direction given by the unit vector f *. 

As it is well known, the system of (R)MHD is not-strictly hyperbolic [21J. This 
means that in some cases, two or more eigenvalues can be equal leading to well studied 
cases of degeneracy (see refs. [II [22]. for the relativistic case). In Type I degeneracy, the 
magnetic field is normal to the propagation direction of the wavefront (i.e., (kB k = 0). In 
Type II degeneracy, (kB k 0, but the eigenvalues associated with, at least, one Alfven 
wave and one magnetosonic wave are degenerate. Leaving aside the particular cases 
associated with both degeneracy types, that will be discussed later, the following list 
compiles the roots of the characteristic equation (1201) . A a (= A*), the right eigenvectors, 
r* [jjj, and their corresponding scalar products, Vf Is in the non-degenerate, general case. 

i) A = A nu n := 0. In this case, V* mn is trivially zero. This eigenvalue is spurious and 
is associated with the fact that although the RMHD system (fTOl) consists of eight 
conservation equations, only seven components of the fluxes are non-trivial. Due 
to the antisymmetric character of the induction equation, the flux of (kB k in the 
(^-direction is identically zero. 

ii) A = Ao := (kV k is the eigenvalue associated with the material waves. 
The corresponding eigenvector is rg = (— k, 0*, y, 0*) T , where k and y are 
thermodynamical derivatives defined at the end of the previous Section, and 0* = 0 
(i = 1,2,3). The scalar product is Vq = 0 and, consequently, the characteristic 
field defined by Ao is linearly degenerate. 

|| The expressions of the eigenvectors have been obtained after tedious algebraic manipulations. They 
can be verified by direct substituting in the eigenvalue equation, ( CkA k — A* A°)r* = 0. 

If For the scalar products V* ± and Vf l± , the partial derivatives of the corresponding eigenvalues with 
respect to the primitive variables, V, have been computed by implicit derivation of the characteristic 
equations for A Q± and A/ ± , respectively, i.e., A = 0 and Afi = 0 (see below). 
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A := Sa 2 - B 2 . (21) 

They define the Alfven waves. Since ( k B k ^ 0, then a ^ 0 and the corresponding 
eigenvectors are 

K ± = (0,4,0,rj) T , (22) 

where r l 2 = a\B l + a 2 v l + a 3 (\ r l = Wa~ 1 (r 2 ( k B k — B l C tk r 2 ). The coefficients 
a p ( p = 1,2,3) are such that v k r 2 = 1, ( k r 2 = — Wa, and B k r k = —v k B k W 2 . The 
scalar products are 

ri = (^r) 4 + (^l) 4 °= (orf + w a = o. ( 23 ) 

in agreement with the linearly degenerate character of the Alfven waves, 
iv) The four eigenvalues A/ ± , A s± , are the roots of the fourth-order polynomial in A, 
■A/ 4 , 

J\f 4 := ( b 2 + pha 2 )a 2 G - W; 2 pha A - a 2 s GB 2 , (24) 

associated with the fast and slow magnetosonic wavespeeds, respectively. Since 
( k B k 0, then a / 0 and the corresponding eigenvectors are 

4 = (ri,rl,r 3 ,rl) T , (25) 

(m = /, s), where 


t'i = pW 3 (phafG + a 2 ) - GB 2 / a S j , 


r 2 = W [GBB 1 + phWa 2 (X m± v l - ( 1 )), 
r 3 = rip/p 2 , 

r\ = phW 3 a ( (A m± v* - C)( k B k - B\X m± aW- 1 - G) 


The scalar products are 

-W™ 

where d, the derivative of A /4 with respect to A at A = A m ,± (m = /, s), J\f 4 {X mj ±), 


d = a 2 s G 2 B(( k B k ) - (G - X m± aW~ l )phWW s - 2 a 4 , 


VI = b 2 G - pha 2 , 


(29) 
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n=\p 


dal 


dp 


p dal 
+ - s 


p de 


W; 


t-E 


b 2 ( 3 


2 pha 2 -|- 


«2(5 - 3a 2 s )B 2 


(30) 


It is interesting to note that d can only be zero in degenerate states, since it is only 
in these states where both A /4 (A) = 0 and Af^(X) = 0 are satisfied simultaneously. 
Let us now discuss the conditions under which the remaining factors in Eq. (1271) can 
become zero. Quantity a is non-zero as far as (kB k 7 ^ 0. On the other hand, it can 
be proven by simple algebraic manipulation of Eqs. M(A) = 0 and A /4 (A) = 0 that 
V* — 0 if and only if the corresponding magnetosonic eigenvalue is also an Alfven 
eigenvalue (i.e., Type II degeneracy). Since we are avoiding degenerate states, and 
G is always non-zero, we shall concentrate on the changes of sign of V%, in order to 
analyze the possible loss of convexity associated with the magnetosonic waves. 
Since in the case of zero magnetic held, the purely relativistic result has to be 
recovered, we shall try now to rewrite expression (l30l) in terms of the relativistic 
fundamental derivative 

P da 2 


Q = 1 + 


- a Q 


(31) 


2 a 2 s dp 

derived in ref. [19]. The sought expression is 

V; = -2a 2 s W 2 £(l~R)g M , (32) 

with Qm, the fundamental derivative for relativistic, magnetized fluids, being 


Qm Q + F, 


where 


F := -W~ 4 
2 s 


C q/ a s ~ R 

1 -R 

B 


(33) 


(34) 


b 2 


I 11 the previous expressions, R := , and c 2 := — stands for the square of the 


£a 2 ’ a ‘ £ 

Alfven velocity. Moreover, in deriving expression il32l) 

following relation among thermodynamical derivatives 


from 

d_ 

dp 


we have used the 


d_ 

dp 


+ 


p d 

~ ~de 




It is important to note that R — 1 if and only if the eigenvalue corresponds to an 
Alfven wavespeed (i.e., it satisfies equation M(A) = 0). Since we are not considering 
degeneracies, we conclude that R 7 ^ 1 for magnetosonic waves and, consequently, 
1) the denominator in the second term of Q m is well defined, and 2) = 0 if and 

only if Qm — 0 . 

The price to pay for using primitive (or conserved) variables in our analysis of 
convexity is the loss of covariance and a dependence of the fundamental derivative 
Qm on kinematics through quantity R. For fast and slow magnetosonic fields, let 
us carry out the analysis of the magnetic correction to the purely hydrodynamic 
(relativistic) fundamental derivative (Eq. (I34|L in the comoving frame (CF, 
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tA = <5 q), which we will name Fqf,™ 
calculation leads to 


Fcp,m 




f, s ) henceforth. A simple algebraic 


(35) 


where <? m are the solutions of the quadratic equation in A 2 , A/ 4 ,cf(A) = 0, namely 


c 


2 

m 


2 (( w2 + a s c u ) ± ( O ^ 2 + 





(36) 


With 4 = ^ kB p 2 and Wj 2 := 1 - cv 2 , u 2 = a 2 + c 2 - a 2 c 2 . 

Taking into account that, for non-degenerate states, a 2 , c 2 E (c 2 , c 2 )0, we have that 
FcF,m > 0 (m = /, s). Now, the transformation of A as a scalar ensures that 
F m > 0 (m — f,s) in any reference frame, with important consequences for the 
influence of the magnetic held on the convexity of the system. 


4. Analysis of convexity in degenerate states 

4-1. Type I degeneracy 

This degeneracy appears in states in which ( k B k = 0. Now, the roots of the characteristic 
equation () 20 lh the right eigenvectors, and the corresponding scalar products have the 
following properties: 

i) A = A nu ii := 0. It is again the spurious eigenvalue analyzed in the previous Section 
associated with the null flux component. "P* u ii is trivially zero. 

ii) The eigenvalue A = Ao := QkV k has multiplicity 5. The corresponding eigenvectors 
are of the form Tq = (r^ai-Eh + a 2 a 3 B l + a 4 4 ) T ; where 4 is an arbitrary 
vector orthogonal to £* and B\ and rq, r 3 and a p , (p = 1,2, 3,4) are functions of 
the primitive variables. Since only the derivative d\/dv k (= 4) is different from 
zero, the scalar product is 

P* = 0(a 1 F t + a 2 Ci) = 0. (37) 

Hence, the characteristic fields defined by Ao arc linearly degenerate. 

iii) A f t ± are the solutions of the quadratic equation in A 

(fr 2 + pha 2 - a 2 s (y k B k ) 2 ) G - Wf 2 pha 2 = 0, (38) 

and are associated with the fast magnetosonic wavespeeds. The explicit expression 
of these eigenvalues when ( k = ( 1 , 0 , 0 ) can be found in ref. [23J. 

The corresponding eigenvectors can be obtained from those of the fast magnetosonic 
eigenvalues in the general case (see Eq. fl25|l ) making ( k B k = 0, i.e., B = a(v k B k ). 

+ In the CF it can be easily proven that AT^cf (o s ) < 0 and A/ 4 j cf(c q ) < 0, implying that both af and 
cf are between the roots of A/ 4 i cf(A) = 0, namely c 2 s , c 2 . 
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P/ ± = 


WjCT 

2 ph 


V* V 


2 5 


where 


S - (v k B k ) 2 

1 - C kV k 


(39) 

(40) 


v; = 


P 


da 2 s 


dp 


+ 


da 2 s 


P 


p de 


W s 2 ({v k B k f - S) 


- 6 2 (3 - a'l) - 2 pha'l + a 2 (5 - 3 a 2 s )(v k B k ) 2 


(41) 


From Eq. (ITS]) , b 2 — (vkB k )' 2 > 0 and then V* is always positive. Hence the possible 
changes of sign of Vf ± coincide with those of V%. Let us note that the expression 
for coincides with that of the general case (Eq. 1301) making B = a{ykB k ). Then, 
proceeding in exactly the same way as in the general case we conclude that the 
fundamental derivative for relativistic, magnetized fluids for Type I degenerate 
states is 


Q 


M,deg I 


S + 2 W - 


-4 

s 


f C d/ a s T?d C gi \ 

\ 1 — Rdegl /’ 


(42) 


where now, i?d eg i = - 7 -• As discussed in the non-degenerate case, i?d eg i 7 ^ 1 , 

and the corresponding factor is Fd eg i > 0 . 

The special case when VkB k = 0 is obtained by making i?d eg i = 0 in the 
previous expression. The same result for this case is obtained through a 
purely hydrodynamical approach (see Appendix in ref. W) by building up a 
thermodynamically consistent EOS incorporating the effects of the magnetic field. 


4-2. Type II degeneracy 

Now, CkB k 7 ^ 0 and, at least, one eigenvalue associated with an Alfven wave and 
an eigenvalue associated with a magnetosonic wave are degenerated. Three cases are 
distinguished. In cases 1 (c a > a s ) and 2 (c a < a s ) one fast or slow magnetosonic 
eigenvalue, respectively, and an Alfven eigenvalue are degenerated. In these cases, as 
discussed in the previous Section, the quantity V{ defined in Eq. (129|) is zero for the 
degenerate eigenvalues and, hence, the corresponding characteristic fields are linearly 
degenerate. When c a = a s (case 3), an Alfven eigenvalue is degenerated with a pair 
(slow and fast) of magnetosonic eigenvalues. Now, quantity d defined in Eq. (1281) is 
also 0, and we have an indetermination in Vf n± (Eq. [22J. In this case, we have checked 
that the dot product of the magnetosonic eigenvectors associated with the degenerated 
fields and the gradient of the Alfven eigenvalue is zero, which means that the degenerate 
characteristic field is again linearly degenerate. 

* As in the non-degenerate case, for the scalar products Vf ± , the partial derivatives of the corresponding 
eigenvalues with respect to the primitive variables, V, have been computed by implicit derivation of 
the characteristic equation ([38]) . 
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The purely (relativistic) hydrodynamical limit can be obtained as a particular case of 
the Type I degeneracy, in which besides having ( k B k = 0 and v k B k = 0, we make b 2 = 0. 
Hence, from Eq. (|33l) . and making R = 0 and c a = 0, we have Q M b 2= 0 = Q. 

We now discuss the classical (magnetized) limits for both degenerate and non¬ 
degenerate states. These limits are obtained by expanding all the quantities in the 
definition of the fundamental derivative in powers of 1/c 2 (c is the speed of light) and 
keeping the leading term. On one hand, the relativistic (non-magnetized) fundamental 
derivative is Q = Q+Oil/c 2 ), where Q is the classical (non-magnetized) counterpart [ T9] . 
On the other hand, R = (CkB k ) 2 /(pc^ cl ) + 0(l/c 2 ), where c m;C i (m = f,s ) is 

Cm ,cl = -J= (a 2 ,cl + B 2 /p± sj (a 2 cl + B 2 /p) 2 - 4a 2 c ff k B k ) 2 /p^j , and a S)Cl stands for 

the classical definition of the sound speed. Hence, we get from Eq. 

(CkB k ) 2 /(pc 2 J 


Gm.cI Q + 2 


■"a, cl 


/< 


cl 




(43) 


In the previous expression, c a . c i stand for the classical definition of the Alfven speed, 

VWp- 

It can be shown that, taking ( k = (1,0,0), the resulting expression of <7 m, c i is 
proportional to the non-linearity factor for the non-linear fields of the (classical) MHD 
system obtained in ref. [20j (see their equation (17)). 

For Type I degenerate states, since R = O (1/c 2 ), 


Qm, 


deg I,cl 


= g + 


"a, cl 

2 Li. 


(44) 


proportional to the corresponding result obtained in ref. |2D] (see their table I). 

Finally for Type II degenerate states, the eigenvalues that are degenerated lead 
to characteristic fields which are linearly degenerate, whereas the (hypothetical) non¬ 
degenerate magnetosonic field (subcases 1 and 2 ) is genuinely non-linear and its 
properties in relation with convexity are governed by the fundamental derivative in 
Eq. (ITT]) . with c mc \ = c S)C i (subcase 1), and c mjC i = Cf. c \ (subcase 2). 


6. Summary and conclusions 

In this paper we have analyzed the influence of the magnetic field in the convexity 
properties of the RMHD equations. To this purpose we have used the approach of 
Lax, based on the analysis of the linearly degenerate/genuinely non-linear nature of 
the characteristic fields. Degenerate and non-degenerate states have been discussed 
separately and the non-relativistic, unmagnetized limits are properly recovered. The 
characteristic fields corresponding to the material and Alfven waves are linearly 
degenerate and, then, not affected by the convexity issue. The analysis of the 
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characteristic fields associated with the magnetosonic waves reveals, however, a 
dependence of the convexity condition on the magnetic field. 

The result is expressed in the form of a generalized fundamental derivative, Eq. (133]) . 
written as the sum of two terms. The first one is the generalized fundamental derivative 
in the case of purely hydro dynamical (relativistic) flow already obtained in ref. [19]. 
The second one contains the effects of the magnetic field. The analysis of this term in 
the comoving frame (extendable to any other reference system given the scalar nature 
of the term) shows that it is always positive leading to the remarkable result that the 
presence of a magnetic field in the fluid reduces the domain of thermodynamical states 
for which the EOS is non-convex, as it happens in the non-relativistic MHD limit [20]. 

We speculate with the possibility that our findings can be relevant in the context 
of massive stellar core collapse. Depending mostly on the pre-collapse stellar magnetic 
field and on the gradient of the rotational velocity, dynamically relevant magnetic fields 
may develop after the core bounce (see, e.g., [251 Ml [27]). Should these magnetic fields 
become as large as the existing numerical models point out, then our results indicate 
that the loss of convexity would be rather limited, if existing at all. However, it is still a 
matter of debate what is the actual level of magnetic field saturation due to the action 
of the Magneto Rotational Instability (MRI; see, e.g., [28] [29]), and hence, whether or 
not the MRI-amplified magnetic field may have the sufficient strength as to impede the 
development of non-convex regions in the collapsed core. It is very likely that under the 
most common conditions (namely, non-rotating or slowly rotating cores), the magnetic 
field will not play central dynamical role in the post-collapse evolution, though it may 
set the time scale for supernova explosions (e.g., [30]). In such cases, we foresee that 
there might exist a range of physical conditions in which a non-convex EOS may render 
a convexity loss in the post-collapse core that cannot be compensated by the growth 
of pre-collapse magnetic fields, e.g., in slowly rotating (including non-rotating) massive 
stellar cores. Addressing this issue by means of numerical simulations is beyond the 
scope of the present work, and will be considered elsewhere. 
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Appendix. Jacobian matrices of the RMHD system in quasi-linear form 

Matrices A 0 and f k A k associated with the system (fTOl) in quasilinear form are: 


( 

w 

CO 

Vo. 

0 

0i 


X°)f 

(A°)S 

(A°)f 

(A°)H 


XX 


(A°)I 

XX B i 

\ 

o i 

o i 

0* 



where 

(A°)f =(1 + e + X )W 2 v\ 

(A°)S = B { Bj + B 2 8 * + hW 2 {8) + 2W 2 v i v j ), 
(A°)f = (p + k)W 2 v\ 

{A°)% = - 8)v k B k - Bi Vj + 2 v i B j , 

{AX =(l + e )W 2 -W + x(W 2 -l), 

(AXj = - B jVk B k + Vj [B 2 + pW 3 {2hW - 1)], 
{AX = P W 2 + k{W 2 - 1), 

{AXbj = - v j v kB k + Bj{2 - l/W 2 ). 


( 

W( k v k 

(C kA k )5 

0 

0 i 

\ 


(CkA k )f 

(CkAXJ 

(C*A fc )f 

XAxy 



(CkA k y p 

(CkA k y vj 

(■ CkA k y e 

(c k A k y Bj 


\ 

o i 

BXj - Sj(kB k 

o i 

8 )C k v k - vXj 

) 


where 

{( k A k )° =pW{W 2 Vj ( k v k + 0), 

(CkA k ) s ; = (1 + e + X )W 2 v\ k v k + x(\ 

{C k A k )5 = {QBj - 8)C l B l )v k B k + B 2 {8]( k v k - (X + v i Q) 

- B\( jVk B k - 2 Vj ( k B k + Bj( k v k ) - v'BjCkB* 

+ phW 2 {8) ( k v k + vXj + 2W 2 v i v j (kV k ), 

{CkA k )f =v i {p + K)W 2 C k v k + CK, 

(C kA k % = CvjV k B k - 8)v k B k Qv l - BX^ 

- v%v k B k + Vj ( k B k - 2Bj( k v k ) - W~ 2 {BXj ~ C% + 8}( k B k ), 
(C kA k y p = (1 + e + x)W 2 CkV k - W(kV k , 

{CkA k ) T vj = -B j C k B k + B 2 C j + pW[C j {hW-l)+v j C k v k W 2 {2hW-l)], 
XkA k ) T e ={p + n)W 2 C k v k , 

(C kA k ) T Bj = 2B 3 Q k v k - vjC k B k - ( jVk B k . 

All the quantities appearing in the definition of the matrices are defined in the body 
of the paper and 0 * = ( 0 , 0 , 0 ) T , 0 j = ( 0 , 0 , 0 ) and 0 * is the null 3x3 matrix. 
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